%load_ext autoreload
%autoreload 2
import warnings
from functools import partial
from optimization import jg_opt
from pisa.administrative_area import AdministrativeArea
from pisa.facilities import Facilities
from pisa.population import WorldpopPopulation
from pisa.population_served_by_isopolygons import get_population_served_by_isopolygons
from pisa.visualisation import (
plot_facilities,
plot_isochrones,
plot_population,
plot_population_heatmap,
)
INFO:matplotlib.font_manager:Failed to extract font properties from /usr/share/fonts/truetype/noto/NotoColorEmoji.ttf: Can not load face (unknown file format; error code 0x2)
INFO:matplotlib.font_manager:generated new fontManager
import os
from dotenv import load_dotenv
# load environment variables from an `.env` file in the local root directory
load_dotenv()
CBC_SOLVER_PATH = os.getenv(
"CBC_SOLVER_PATH"
) # path to the cbc executable (e.g. /opt/homebrew/bin/cbc)
Pisa Showcase, Example Notebook
Define Administrative Area
Let op: the administrative area is a country because of the package that we use for fetching the data (gadm).
The naming is confusing: we want the administrative area to be a subset of a country (in this case, the administrative area should be Baucau).
Should be easy to fix once we change gadm to pygadm (see issue #59 on GitHub)
timor_leste = AdministrativeArea(country_name="Timor-Leste", admin_level=1)
# these are the boundaries of Baucau
# type: Polygon
baucau = timor_leste.get_admin_area_boundaries("Baucau")
INFO:pisa.administrative_area:Validating country name: Timor-Leste
INFO:pisa.administrative_area:Country name 'Timor-Leste' validated successfully
INFO:pisa.administrative_area:Retrieving boundaries of all administrative areas of level 1 for country Timor-Leste
0%| | 0.00/6.25M [00:00<?, ?B/s]
1%| | 39.9k/6.25M [00:00<00:18, 334kB/s]
3%|▎ | 209k/6.25M [00:00<00:06, 962kB/s]
14%|█▍ | 889k/6.25M [00:00<00:01, 3.24MB/s]
29%|██▉ | 1.83M/6.25M [00:00<00:00, 5.33MB/s]
38%|███▊ | 2.38M/6.25M [00:00<00:00, 4.84MB/s]
46%|████▌ | 2.88M/6.25M [00:00<00:00, 4.50MB/s]
54%|█████▍ | 3.40M/6.25M [00:00<00:00, 4.62MB/s]
65%|██████▍ | 4.06M/6.25M [00:00<00:00, 4.72MB/s]
73%|███████▎ | 4.58M/6.25M [00:01<00:00, 4.81MB/s]
84%|████████▍ | 5.24M/6.25M [00:01<00:00, 5.26MB/s]
94%|█████████▍| 5.90M/6.25M [00:01<00:00, 5.32MB/s]
100%|██████████| 6.25M/6.25M [00:01<00:00, 4.65MB/s]
Facilities
Get existing facilities (in our case, hospitals) from OSM
# Suppress user warning about geometry in geographic CRS. Centroid is calculated over a single facility (e.g. a hospital),
# so distances are very small and projection isn't necessary
warnings.filterwarnings(
"ignore",
message="Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect",
category=UserWarning,
)
hospitals_df = Facilities(admin_area_boundaries=baucau).get_existing_facilities()
INFO:pisa.facilities:Retrieving existing facilities with tags {'amenity': 'hospital', 'building': 'hospital'} using OSM.
INFO:pisa.facilities:Successfully retrieved existing facilities from OSM.
plot_facilities(hospitals_df, baucau)
Make this Notebook Trusted to load map: File -> Trust Notebook
Estimate potential locations for new facilities
potential_hospitals_df = Facilities(
admin_area_boundaries=baucau
).estimate_potential_facilities(spacing=0.05)
Get population
In this example, from WorldPop
# TODO: why "ID" column? extra index for what reason? Maybe becomes apparent later?
population_gdf = WorldpopPopulation(
admin_area_boundaries=baucau, iso3_country_code=timor_leste.get_iso3_country_code()
).get_population_gdf()
population_gdf.head()
| ID | longitude | latitude | population | geometry | |
|---|---|---|---|---|---|
| 0 | 0 | 126.14917 | -8.67333 | 0.10 | POINT (126.14917 -8.67333) |
| 1 | 1 | 126.14917 | -8.67250 | 0.09 | POINT (126.14917 -8.6725) |
| 2 | 2 | 126.15000 | -8.67333 | 0.10 | POINT (126.15 -8.67333) |
| 3 | 3 | 126.15000 | -8.67250 | 0.09 | POINT (126.15 -8.6725) |
| 4 | 4 | 126.15000 | -8.67167 | 0.09 | POINT (126.15 -8.67167) |
plot_population_heatmap(population_gdf, baucau)
Make this Notebook Trusted to load map: File -> Trust Notebook
plot_population(population_gdf, baucau, random_sample_n=1000)
Make this Notebook Trusted to load map: File -> Trust Notebook
Calculate isopolygons
Here we make some choices:
distance type
distance values
mode of transport
Valid values for constants are in the script pisa.constants
DISTANCE_TYPE = "length"
DISTANCE_VALUES = [2000, 5000, 10000]
MODE_OF_TRANSPORT = "driving"
Using OSM
For this, we need to get the road network from osmnx
from pisa.osm_road_network import OsmRoadNetwork
road_network = OsmRoadNetwork(
admin_area_boundaries=baucau,
mode_of_transport=MODE_OF_TRANSPORT,
distance_type=DISTANCE_TYPE,
).get_osm_road_network()
INFO:pisa.osm_road_network:OSM road network set with parameters network_type 'drive' and distance_type 'length'
Calculate isopolygons for existing facilities
from pisa.isopolygons import OsmIsopolygonCalculator
isopolygons_existing_facilities = OsmIsopolygonCalculator(
facilities_df=hospitals_df,
distance_type=DISTANCE_TYPE,
distance_values=DISTANCE_VALUES,
road_network=road_network,
).calculate_isopolygons()
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
isopolygons_existing_facilities.head()
| ID_2000 | ID_5000 | ID_10000 | |
|---|---|---|---|
| osmid | |||
| 1777614876 | POLYGON ((126.59878421714033 -8.54764218472667... | POLYGON ((126.59761276961231 -8.57116537023302... | POLYGON ((126.56068305448515 -8.56751807760911... |
| 1777614896 | POLYGON ((126.65907104033573 -8.62868878467725... | POLYGON ((126.6311023795325 -8.646327183432366... | POLYGON ((126.63106042126435 -8.64641589673682... |
| 1777615024 | POLYGON ((126.59499222126435 -8.47806089673682... | POLYGON ((126.55950911149382 -8.48560952854904... | POLYGON ((126.57412141532275 -8.51524054033573... |
| 1777615035 | POLYGON ((126.68354761081066 -8.63974323564052... | POLYGON ((126.68357328721343 -8.63970069138097... | POLYGON ((126.6311023795325 -8.646327183432366... |
| 1777615036 | POLYGON ((126.6803646 -8.578027, 126.680359784... | POLYGON ((126.6650524795325 -8.602129783432366... | POLYGON ((126.65894556961231 -8.62895407023302... |
plot_isochrones(isopolygons_existing_facilities, baucau)
Make this Notebook Trusted to load map: File -> Trust Notebook
Calculate isopolygons for potential facilities
isopolygons_potential_facilities = OsmIsopolygonCalculator(
facilities_df=potential_hospitals_df,
distance_type=DISTANCE_TYPE,
distance_values=DISTANCE_VALUES,
road_network=road_network,
).calculate_isopolygons()
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:234: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:239: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.
cylinders = edges_gdf.buffer(edge_buffer).geometry
Prepare optimization data
Preparing the variables that go into the cell below
population_served_current = get_population_served_by_isopolygons(
grouped_population=population_gdf, isopolygons=isopolygons_existing_facilities
)
current = {DISTANCE_TYPE: population_served_current}
population_served_potential = get_population_served_by_isopolygons(
grouped_population=population_gdf, isopolygons=isopolygons_potential_facilities
)
potential = {DISTANCE_TYPE: population_served_potential}
population_count = population_gdf.population.values
BUDGET = [
5,
20,
50,
] # budget for the optimization in terms of how many locations can be built
cbc_optimize = partial(jg_opt.OpenOptimize, solver_path=CBC_SOLVER_PATH)
values, solutions = jg_opt.Solve(
household=population_count,
current=current,
potential=potential,
accessibility=DISTANCE_TYPE,
budgets=BUDGET,
optimize=cbc_optimize,
type="ID",
)
WARNING:pyomo.opt:Failed to create solver with name '_cbc_shell':
Failed to set executable for solver cbc. File with name=D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe either does not exist or it is not executable. To skip this validation, call set_executable with validate=False.
Traceback (most recent call last):
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 148, in __call__
opt = self._cls[_name](**kwds)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/solvers/plugins/solvers/CBCplugin.py", line 88, in __init__
super(CBCSHELL, self).__init__(**kwds)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/solver/shellcmd.py", line 66, in __init__
self.set_executable(name=executable, validate=validate)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/solver/shellcmd.py", line 115, in set_executable
raise ValueError(
ValueError: Failed to set executable for solver cbc. File with name=D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe either does not exist or it is not executable. To skip this validation, call set_executable with validate=False.
WARNING:pyomo.opt:Failed to create solver with name 'cbc':
Attempting to use an unavailable solver.
The SolverFactory was unable to create the solver "_cbc_shell"
and returned an UnknownSolver object. This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "set_problem_format").
The original solver was created with the following parameters:
executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
type: _cbc_shell
_args: ()
options: {}
Traceback (most recent call last):
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 148, in __call__
opt = self._cls[_name](**kwds)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/solvers/plugins/solvers/CBCplugin.py", line 50, in __new__
opt.set_problem_format(ProblemFormat.cpxlp)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 113, in __getattr__
self._solver_error(attr)
File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 116, in _solver_error
raise RuntimeError(
RuntimeError: Attempting to use an unavailable solver.
The SolverFactory was unable to create the solver "_cbc_shell"
and returned an UnknownSolver object. This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "set_problem_format").
The original solver was created with the following parameters:
executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
type: _cbc_shell
_args: ()
options: {}
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call last)
Cell In[20], line 9
1 BUDGET = [
2 5,
3 20,
4 50,
5 ] # budget for the optimization in terms of how many locations can be built
7 cbc_optimize = partial(jg_opt.OpenOptimize, solver_path=CBC_SOLVER_PATH)
----> 9 values, solutions = jg_opt.Solve(
10 household=population_count,
11 current=current,
12 potential=potential,
13 accessibility=DISTANCE_TYPE,
14 budgets=BUDGET,
15 optimize=cbc_optimize,
16 type="ID",
17 )
File ~/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/optimization/jg_opt.py:362, in Solve(household, current, potential, accessibility, budgets, optimize, type)
360 I = np.unique(list(IJ.keys()))
361 J = np.unique(np.concatenate(list(IJ.values())))
--> 362 optimization = optimize(
363 household,
364 I,
365 J,
366 IJ,
367 [max(budgets)],
368 parsimonious=True,
369 maxTimeInSeconds=5,
370 mipGap=1e-15,
371 )
372 optimization["nof"] = [len(s) for s in optimization.solution]
373 coverage = (optimization.value / household.sum() + percent_covered).to_frame()
File ~/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/optimization/jg_opt.py:126, in OpenOptimize(w, I, J, IJ, budget_list, parsimonious, maxTimeInSeconds, mipGap, trace, solver, solver_path)
124 result.at[p, "modeling"] = pc() - start
125 start = pc()
--> 126 solver_result = solver.solve(M, tee=trace)
127 result.at[p, "solving"] = pc() - start
128 result.at[p, "value"] = int(
129 np.ceil(M.weighted_coverage() - np.finfo(np.float16).eps)
130 )
File ~/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py:99, in UnknownSolver.solve(self, *args, **kwds)
97 def solve(self, *args, **kwds):
98 """Perform optimization and return an SolverResults object."""
---> 99 self._solver_error('solve')
File ~/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py:116, in UnknownSolver._solver_error(self, method_name)
115 def _solver_error(self, method_name):
--> 116 raise RuntimeError(
117 """Attempting to use an unavailable solver.
118
119 The SolverFactory was unable to create the solver "%s"
120 and returned an UnknownSolver object. This error is raised at the point
121 where the UnknownSolver object was used as if it were valid (by calling
122 method "%s").
123
124 The original solver was created with the following parameters:
125 \t"""
126 % (self.type, method_name)
127 + "\n\t".join("%s: %s" % i for i in sorted(self._kwds.items()))
128 + "\n\t_args: %s" % (self._args,)
129 + "\n\toptions: %s" % (self.options,)
130 )
RuntimeError: Attempting to use an unavailable solver.
The SolverFactory was unable to create the solver "cbc"
and returned an UnknownSolver object. This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "solve").
The original solver was created with the following parameters:
executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
type: cbc
_args: ()
options: {'threads': 8, 'sec': 5.0, 'allowableGap': 1e-15}
values
| 10000 | 5000 | 2000 | |
|---|---|---|---|
| 5 | 0.385968 | 0.256445 | 0.126371 |
| 20 | 0.393635 | 0.266618 | 0.133419 |
| 50 | 0.393635 | 0.268916 | 0.134770 |
solutions
| 10000 | 5000 | 2000 | |
|---|---|---|---|
| 5 | [27, 43, 16, 10, 17] | [27, 22, 25, 9, 43] | [39, 43, 25, 7, 22] |
| 20 | [27, 43, 16, 10, 17, 22, 18, 35, 37, 14, 38, 0... | [27, 22, 25, 9, 43, 16, 18, 42, 50, 38, 39, 28... | [39, 43, 25, 7, 22, 27, 24, 26, 13, 32, 10, 16... |
| 50 | [27, 43, 16, 10, 17, 22, 18, 35, 37, 14, 38, 0... | [27, 22, 25, 9, 43, 16, 18, 42, 50, 38, 39, 28... | [39, 43, 25, 7, 22, 27, 24, 26, 13, 32, 10, 16... |